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DETERMINATION  OF  AHRQDTNAM1C  COEFFICIENTS  USINQ  ACCELEROMETER  RECORDS 

FROM  A PUNE  IAWINO  BOMB 


ABSTRACT 

The  presently  used  Qavre  drag  functions  are  not  well  suited  to 
the  streamlined  bombs  dropped  from  high  speed,  high  altitude  airplanes. 
One  method  of  obtaining  better  drag  functions  is  from  accelerometers 
mounted  inside  the  bombs . This  report  develops  a method  of  reducing 
such  accelerometer  data  and  applying  existing  exterior  ballistics 
theories  to  obtain  not  only  K^,  but  K^,  and  spin  as  well. 

Five  bombs  equipped  with  accelerometers  were  dropped  at  White  Sands 
Proving  Ground,  New  Mexico  during  February  and  March  1952.  Although 
the  results  were  not  too  satisfactory  they  indicate  the  method  is 
accurate  aad  practical  if  the  following  can  be  obtained:  (a)  an 

electrical  system  giving  a smooth  record  of  accelerations  with  known 
limits  of  error,  (b)  accelerometers  of  two  different  ranges  for  small 
relative  errors  in  drag,  (c)  accurate  meteorological  data  near  the  time 
of  drop. 
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INTRODUCTION 


In  recent  years  higher  airplane  speeds  and  altitudes  hare  resulted 
in  bomb  speeds  greater  than  that  of  sound  over  an  appreciable  part  of 
the  trajectory.  To  ensure  good  flight  at  such  speeds,  bombs  have  been 
made  more  streamlined,  with  longer  tails  and  fins  chosen  for  aerodyna- 
mic qualtities.  The  drag  coefficients  of  such  bombs  are  not  well  repre- 
sented by  the  coefficient  corresponding  to  the  Qavre  drag  function)  the 
new  bombs  have  less  drag  in  the  subsonic  region  and  a steeper  rise  at 
the  critical  Mach  number,  which  gives  larger  drags.  The  result  is  that 
the  Bomb  Ballistic  Reduction  Tables,  based  on  the  Qavre  drag  function, 
can  only  be  applied  to  the  new  bombs  by  violent  means.  It  is  necessary 
to  have  different  ballistic  coefficients  for  range  and  time  of  flight, 
and  to  vary  them  with  launching  conditions. 

The  fast  electronic  computing  machines  now  available  make  the 
computation  of  a bombing  table  from  a given  drag  function  rather  simple. 

It  then  appears  (lesirable  to  compute  the  table  for  each  bomb  directly, 
using  its  own  drag  function. 

Three  methods  of  getting  the  necessary  aerodynamic  information  are 
available!  Wind  tunnel  tests,  spark  range  firings,  and  drops  of  bombs. 
Unfortunately,  drag  measurements  at  supersonic  speeds  are  hard  to  make 
in  a tunnel,  and  of  lower  accuracy  than  seems  desirable.  Also,  wind 
tunnel  and  spark  range  tests  must  have  Reynolds  numbers  far  from  those 
of  actual  flight,  because  small  models  must  be  used)  and  the  effect  of 
different  Reynolds  numbers  ("scale  effect")  is  not  well  understood.  Drops 
of  full-scale  models,  properly  instrumented,  seem  to  offer  these  advantages 
(1)  data  are  obtained  from  duplicates  of  the  actual  bomb,  (even  the 
surface  finish,  for  example,  is  the  same)  dropped  under  tactical  con- 
ditions) (2)  measurements  of  drag,  lift,  moment  and  damping  rate  can  be 
obtained  from  the  same  drop)  (3)  data  for  the  construction  of  a bombing 
table  should  be  obtainable  from  a small  number  of  drops,  without  ex- 
tensive range  bombings)  (U)  the  behaviour  of  the  actual  bomb  under  all 
expected  conditions  of  flight  can  be  examined  directly,  which  gives  a 
satisfaction  that  calculations  from  models  do  not  provide. 

It  was  accordingly  decided  to  drop  demolition  bombs  of  the  new 
family,  T-5h,  T-55,  and  T-J>6,  at  White  Sands  Proving  Qround. 

A sketch  of  the  bomb  showing  approximate  locations  and  positions 
of  the  accelerometers  is  shown  in  Pig.  10. 
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Description  of  the  apparatus.  At  the  center  of  gravity^-  of  each  boob 
were  placed  three  accelerometers*  The  first  vas  placed  along  the  axle 
of  symmetry,  and  indicated  axial  drag  accelerations  from  sero  to  2g. 

Its  output  vent  on  Channel  1 of  the  telemetering  apparatus.  The  second 
and  third  accelerometers  were  on  the  axis-'-  at  right  angles  to  the  first 
and  to  each  other,  to  measure  accelerations  normal  to  the  axis.  Their 
outputs  vent  on  Channels  2 and  3*  Approximately  eighteen  inches  behind 
the  c.g.,  and  on  the  axis,  were  two  more  accelerometers  on  Channels  U 
and  £,  parallel  respectively  to  the  second  and  third.  As  far  back  as 
they  could  be  conveniently  located  in  the  tail  (h  ft.  from  the  c.g.  in 
the  text)  vere  tvo  more  accelerometers  on  the  axis,  respectively  parallel 
to  the  second  and  third,  using  Channels  6 and  7»  The  fourth  and  fifth 
accelerometers  vere  carried  as  insurance.  Estimates  indicated  that  if 
undesirably  large  yaws  developed  near  the  speed  of  sound,  the  rearward 
accelerometers  might  be  overloaded.  The  fourth  and  fifth  vere  put  in  to 
measure  such  yaws,  but  none  vas  observed.  All  accelerometers  but  the 
first  had  the  range  -lg  to  +lg.  The  first  five  accelerometers,  the 
telemetering  apparatus,  the  power  supplies,  and  a Sphered op  apparatus 
for  position-time  data  were  mounted  on  a steel  tray  which  could  be  slid 
into  the  inert  loaded  boob  before  the  tail  vas  mounted. 

Ground  equipment  Included  ballistic  cameras  to  get  position  and 
velocity  at  launch;  Askania  theodolites  far  position-time  data  during 
flight,  in  case  the  Sphered op  didn't  work  veil;  Bowen-Knapp  cameras 
to  cover  the  last  thousand  feet  of  the  trajectory,  to  give  the  Sjphere- 
dop  data  a veil-determined  origin;  telemetry  receiving  stations;  radar 
tracking  equipment;  and  weather  balloons. 

The  telemetering  and  Sphered op  apparatus  vere  designed,  installed, 
and  tested  by  the  Ballistic  Measurements  Laboratory  of  BRL;  the  ground 
equipment  vas  operated  by  the  Flight  Determination  Laboratory  of  WSPG; 
the  meteorological  data  vere  furnished  by  the  AAF  weather  station  at 
WSPG.  The  airplane  and  crew  were  assigned  by  the  Aberdeen  Boobing 
Mission  at  Edwards  AFB  from  the  machines  and  crews  stationed  there  by  • 
the  AAF  for  boobing  tests.  The  telemetering  records  vere  read  by  the 
Flight  Determination  Laboratory,  and  the  Sphered  op  data  vere  reduced  by 
the  Ballistic  Measurements  Laboratory. 

Accelerometer  records.  Each  accelerometer  put  out  a voltage  be- 
tveen  0 and  volts,  linear  in  the  acceleration  affecting  it.  For  drag, 
the  output  vas  sero  for  zero  acceleration;  for  the  other  accelerometers 
the  output  vas  2.5  volts  at  zero.  The  output  controlled  the  pulse  width 
in  a pulse  width  frequency  modulation  telemetering  circuit.  Bach 
accelerometer  reading  was  sampled  about  20  times  per  second.  (The 
speed  of  the  coonutator  varied  somewhat.)  The  ground  station  shoved  the 


1.  The  actual  size  of  the  accelerometers,  and  the  small  variations  of 

position  of  c.g.  from  bomb  tft  bomb; tfctoe  Sts s only  applfximat ely 
correct.  No  appreciable  error  vas  introduced  by  this  approximation. 
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output  of  each  acoelerometer  u the  length  of  a line  an  an  oscilloscope. 
The  oscilloscope  vas  photographed  by  a moving  picture  camera.  Since  the 
telemetering  apparatus  had  thirteen  channels,  six  were  used  to  send 
reference  roltagss  of  0,  1,  2,  3,  1*,  $ volts.  (A  temperature  indicator 
vas  used  instead  of  3 volts  on  some  drops.) 

The  film  vas  read  by  a Homograph,  vhich  measured  the  length  of  each 
line  by  means  of  a photoelectric  cell  and  marked  a corresponding  point 
on  paper.  The  Homograph  adjusted  itself  automatically  to  the  zero 
reference  voltage  line,  and  manual  control  vas  used  to  try  to  fit  another 
reference  voltage,  k or  £ volts. 

A section  of  film  vas  read  at  BRL,  using  an  ordinary  reader  and 
interpolating  between  the  nearest  two  reference  voltages,  instead  of  be- 
tween 0 and  h or  5 volts.  The  resulting  record  was  smoothed;  but  In 
view  of  the  methods  to  be  used  in  determining  the  aerodynamic  coeffic- 
ients, the  labor  of  reading  all  the  film  did  not  seem  worth  while. 

Conduct  of  the  tests.  Because  of  various  delays,  conflict  with  other 
programs  arose,  and  only  five  bombs  were  dropped.  The  available  air- 
plane, a B-29,  could  not  attain  the  desired  speed  and  altitude.  The 
bomb  vas  slung  below  the  aircraft,  from  which  the  boob  bay  doors  had 
been  removed,  so  that  most  of  the  boob  vas  outside  the  bay.  (Two 
3000  lb.  T -55  bombs  were  carried  at  once.)  The  slings  were  designed  to 
release  the  bomb  with  practically  no  disturbance  in  yaw  or  spin.  It 
was  not  possible  to  cock  the  nose  up,  as  this  would  have  brought  the 
tail  too  close  to  the  runway.  For  an  initial  yaw,  the  angle  of  attack 
of  the  airplane  and  the  curved  airstream  near  the  fuselage  were  depended 
on.  The  accelerometers  normal  to  the  axis  pointed  down  at  about  U5* 
from  the  vertical.  The  instruments  in  the  boob  were  connected  to  the 
airplane 's  power  supply  through  a pull-out  plug. 

Since  the  airfield  at  WSPG  was  small,  it  vas  necessary  to  fly  from 
Edvards  Air  Force  Base,  After  take-off,  the  airplane  flew  over  the 
Instrument  building  and  the  accelerometers  (except  the  first)  vers  read 
to  be  sure  they  indicated  about  lg  cos  l£°.  The  airplane  then  flew  to 
WSPG,  reporting  its  arrival  over  Albucperque . (From  Albuquerque  to  WSPG 
the  course  vas  approximately  the  desired  bombing  course.)  As  the  air- 
plane approached  WSPG  its  course  vas  plotted  on  the  radar  plotting 
board  and  it  vas  talked  on  to  the  desired  line  of  flight.  The  bombardier 
picked  up  his  target,  a specially  prepared  circle  200  ft.  in  diameter. 

The  airplane  than  made  a large  circle,  giving  time  far  the  instruments 
in  the  bomb  to  be  checked  by  the  ground  stations.  About  two  minutes 
before  the  drop  the  instruments  vers  connected  to  the  bomb's  internal 
power  supply  and  their  operation  checked  again.  Last  minute  corrections 
to  the  lias  of  flight  were  made  as  the  airplane  approached  its  straight 
bombing  run,  and  the  release  of  the  bcmb  was  controlled  by  the  bombardier 
with  a standard  optical  sight.1 2 


1.  As  miggssted  by  S.  S.  Martin,  of  BRL. 

2.  After  Drop  No.  3,  of  a 3000  lb.  bomb,  one  engine  failed,  and  the 
second  bomb  of  the  pair  was  dropped  with  only  three  engines  working. 
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Remarks.  The  reaulta  of  the  experiment  aa  shown  in  Section  U ware 
disappointing.  Nevertheless,  these  tests  did  hare  a amber  of  useful 
results.  First,  it  was  shown  that  the  T-55  aid  T-56  bsrnbs  flew  well, 
developing  no  objectionable  yaws  at  speeds  near  that  of  sound*  Second, 
the  Spheredop  apparatus  gave  good  position  tine  data.  Third,  aetheds 
of  analysing  the  data  were  developed  which,  judging  from  the  results  of 
these  drops,  will  give  cpite  accurate  values  of  the  aerodynanic  coeffi- 
cients if  smooth  data  can  be  obtained. 


1.  METHOD  AMD  THEORY 

In  this  section,  we  derive  the  equations  from  which  several  of  the 
aerodynamic  coefficients  and  the  spin  can  be  determined.  The  coeffi- 
cients corresponding  to  the  drag  D,  lift  L,  normal  force  N,  restoring 
moment  M,  and  damping  moment  H will  be  defined  by  the  following 
equations: 

(1.1)  D - KjjpdV, 

(1.2)  L - , 

(1.3)  N - KjjpdVfi  , 

(1.U)  M - -KjjpdVs  , 

(1.5)  H - -Kjjpd^uo, 

where  p is  the  density  of  the  air  in  lbs  ./ft  ? , d is  the  diameter  of  the 
bomb  in  ft.,  u is  the  air  speed  of  the  bomb  in  ft./sec.,  6 is  the  angle 
of  yaw  in  radians,  and  oo  the  angular  velocity  of  the  longitudinal  axis 
of  the  bomb  in  rad  ./sec. 

Because  the  bombs  had  uncanted  fins  and  were  suspended  nearly 
horizontally,  being  released  with  negligible  angular  velocity,  we  shall 
make 

Assumption  1.  The  spin  rate  is  mall. 

Assumption  2.  The  yaw  is  small  and  nearly  planar. 

The  first  of  these  is  borne  out  by  the  data  and  the  second  is  a conse- 
-qcu.rce  of  thw  first.  Furthermore,  assumption  2 justifies  making 

Assumption  3.  The  axial  drag  D.  represents  the  total  drag  within  the 
accuracy  of the  experiment . * 

This  assumption  in  borne  out  by  calculations  from  the  data. 
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How  Ba  - ma^  which  together  with  (1*1)  and  assumption  3 give 

(1.6)  Kjj  - au^/pdV 

where  n is  the  mass  of  the  projectile  in  pounds  and  a.  is  the  accelera- 
tion as  recorded  on  channel  j,  J ■ 1,  2,  ...,  7.  J 

In  Section  3 it  will  be  shown  using  the  results  of  McShane  [lj 
under  assumptions  1,  2,  7,  8 that  the  differential  equation  of  yaw  6 is 


(1.7) 

6"  ♦ 2(a  + u'/2u)6»  + p26  1 

and  the  yaw  is 

given  closely  by 

(1.8) 

8 - f/u^u  e”a*(cei^8  + ce 

where 

(1.9) 

a - (pd2/^2)^  ♦ k\) 

(1.10) 

P • <V^dK^/mk2 

and  primes  indicate  derivatives  with  respect  to  arclength  s along  the 
trajectory,  u - An  . - 0,  c aod  Z are  copier  conjagat^tu.ta, 

i • V-l,  B is  the  transverse  moment  of  inertia,  and  k ■ "y  B/md  , the 
radius  of  gyration  in  calibers. 

Let  62  denote  the  two  complex  terms  of  6 in  (1.8)  and  <p^ 
the  arguments  of  6^,  6g  respectively.  Then 


(1.11) 

6 " ®1  + ®2  * 

Cf  1 ",  P«  + T , 

(1.12) 

1^2  " “P®  " T » 

where  y is  the  argaent  of  c.  The  yaw  is  shown  below  in  the  complex 
plane. 


• f 


From  equations  (1.12)4 

d.U)  fi  - 92  ■ tfv  ft  - -fr 

6 is  aero  whenever  ■ n/2  + up  for  p an  integer.  Let  X/2  be  the  arc 
length  between  successive  zero! of  6 and  call  X the  wave  length.  Using  (1.12) 

(P(s  + X/2)  * y)  - (P»  + r)  ■ * 
or 

(1.14)  pX  - 2n. 

Substituting  (l.lU)  into  (1.10)  yields 

(1.15)  ^ - yVmlcZAZpd 

We  shall  call  the  points  midway  in  arclength  between  zeros  of  6 
the  "midarcpointa"  between  zeros  of  6.  The  points  where  <p^  » 0, 

n,  2n,  ...  are  such  points.  At  the  midarcpointa  between  zeros  of  6, 
we  have  62  ■ 6^, 

In  tents  of  tine  derivatives  (1.7)  becomes 

(1.16)  6 ♦ 2ol  + p2u26  • 0. 

Trm  (1.8),  (l.U)  we  obtain 

(1.17)  i - -(«  ♦ u«/2u)uC  + ipu^  - 62) 
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•ad  »t  the  midarcpoints  between  seros  «f  6, 


ao£  - -2a(«  * u'/2u)u6 

2 2 

which  we  coapare  with  jru  6 in  (1,16).  .For  the  T£j>  end  T$6  bombs  with 
airspeed  £00  ft  ./sec.  we  obtained 

2a6  ~10*U6  , 

pVfi  ~ 2.56  , 

6 -5&1  . 

In  general  we  can  make  the 

Aeaumptien  U«  At  the  nldarcpoints  between  seres  of  6,  6 can  be  neglected 

In  (1.16). 

At  the  seros  ef  6,  62  ■ -6^  and  (1.16),  (1.17)  yield 

• • « 

6 - -2o6  , 


6 - 2i0u6^  , 


or 


6 » -UioPu61  . 

Fair  the  T55  and  T56  boobs  with  airspeed  500  ft. /sec.  we  obtained1 
6 ~ -2i  x 10"^  . 


In  general  we  can  make  the 

## 

Aggtsaption  5.  The  seros  of  6 and  6 coincide  within  the  error  measure- 
meat. 

Using  assumptions  U and  5 we  can  make  the 

^•sumption  6.  At,  the  midarcpoints  of  the  seros  of  6,  6 can  be  neglected 
Using  assumption  6 we  hare 

(1.18)  S - -u2p26 
at  the  midarcpoints  between  seres  ef  6. 

1 For  the  ralues  ef  the  ingredients  in  this  estimate  see  the  discussion 
of  Assumption  8 in  Section  3, 
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Define  the  resultist  r of  the  transverse  accelerations  (at  the  rear- 
ward acceleraaetera)  due  to  yawing  motion  by 

(1.19)  r - Va61 2  ♦ a?2  - . 

Because  accelerations  a2  and  a,,  act  through  the  center  of  gravity  the 

second  term  on  the  right  in  (1.19)  eliminates  the  effect  of  lift  upon 
the  transverse  accelerations.  If  accelerometers  6 and  7 are  at  a dis- 
tance^- h from  the  center  of  gravity, 

(1.20)  r - h ( 6 { . 

Equations  (1.8),  (1.18),  (1.20)  combine  to  give 

r - h|^VV/V«  | ♦ &-1H 

Let  C be  such  a positive  number  that 

(1.21)  r - Cu3/V“ 

1 

• • 

for  the  midarcpoints  between  the  zeros  of  6 or  between  the  minima  of  r. 
Equation  (1.21)  will  be  used  to  determine  a. 

Now  |N|  ■ + a^2  which  together  with  (1.3)  yields 

(1 .22  ) Kjj  ■ m Va^+~a^/pd2u2 1 6 1 

Assumption  1 together  with  (l.l),  (1.2)  and  (1.3)  yield  the  usual 
formula 

(1.23)  ^ - Kjj  - Kjj  . 

We  can  new  determine  by  solving  (1.9) 

(1.210  Kh  - (-1^  + 2ma/pd2)k2 

This  analysis  should  be  sligitly  modified  if  the  bomb  spins  slowly, 
but  the  position  and  size  of  the  maxima  and  minima  of  r will  not  be 
affected  within  the  errors  of  measurement. 


1 Accelerometers  6 and  7 may  be  replaced  by  it  and  in  (1.19)  and  in 

the  definition  of  h. 
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To  determine  the  axial  apin  choose  non-rotating  reference  axes  so 
that  the  transverse  acceleration  (at  accelerometers  6,  7)  is  parallel  to 
the  axis  of  the  abcissa.  Let  a be  the  signed  magnitude  of  this  trans- 
verse acceleration1* 


Let  6,  7 be  such  orthogonal  axes  that  a^,  a^  are  the  components  of  the 

above  acceleration  in  directions  6,  7 respectively.  Let  y be  the  angle 
from  the  positive  abcissa  axis  to  axis  6.  Then 


(1.25) 
and 

(1.26) 

Let 

(1.27) 


a^  » a cosy  , 
a^  ■ -a  siny 


2 


a 


ViVs2  * */ 
-7/I V*62  * *1  . 


where  the  sign  is  chosen  as  that  of  a.  For  the  latter  purpose  we  must 

2 

keep  track  of  the  minima  of  a (they  were  practically  zero)  which  should 
be  the  same  as  the  seros  of  S. 


1 The  year  was  practically  plane  in  the  part  of  the  trajectory  con- 

sidered. 
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Combining  (1.25)  and  (1.26)  i 

(ll28) 

( cosy  - 
siny  ■ 

or 

* 

( y ■ arc 

(1.29) 

\ y-  arc 

In  order  to  specify  uniquely 

c« 

an 


we  specify  an  initial  y , require  y be  continuous,  and  at  extrema  of  Cs, 
Sn  use  Sn,  Ca  respectively  to  determine  whether  y is  increasing  or  de- 
creasing. 


The  rate  of  spin  is  given  by  \p 


2.  COMPUTATIONAL  AND  FITTING  PROCEDURES 

Required  data.  To  determine  the  aerodynamic  coefficients  and  spin  of  a 
bomb  using  the  equations  of  Section  1,  it  is  necessary  to  have  avail- 
able the  physical  data  of  the  bomb  (moments  of  inertia,  mass,  caliber), 
meteorological  data  (wind  velocity  and  air  density  and  temperature  vs 
altitude),  trajectory  data  (altitude  and  velocity  vs  time),  and  the 
readings  of  the  accelerometers  vs  time.  The  air  speed  u can  be  obtained 
by  correcting  the  velocity  of  the  bomb  with  respect  to  the  ground  by  the 
wind  velocity.  For  the  trajectory  referred  to  air,  the  arclength,  s, 
as  a function  of  time,  t,  can  be  obtained  by  integrating  the  airspeed  u 
with  respect  to  t. 

Damping  rate  and  aerodynamic  coefficients . The  drag  coefficient  can  be 
obtained  by  substitution  in  (1.6)  • 

Equation  (1.19)  and  the  readings  of  the  accelerometers  can  be  used 
to  determine  the  r vs  t function.  The  r vs  s function  can  then  be 
plotted  and  the  values  of  s for  the  minima  of  r (corresponding  to  zeros 
of  6 and  2$)  can  be  measured.  In  accord  with  assumption  $ the  measured 
arclengths  „Op  -o  2*  •••>  ^ n+1  from  a convenient  origin  to  the  minima 

of  r will  be  used  to  obtain  fitted  values  Sp  S2,  ...,  SQ+1  which  differ 

by  a uniform  interval  X/2 . The  least  square  fit  to  the  measured  inter- 
vals is  given  by 

(2.1)  X/2  - tin+1  -^jVn  . 

Nov  0 and  Kg  can  be  determined  using  (2.1)  in  (1.14)  and  (1.15), 
respectively.  Setting 

(2.2)  8i  " ®l  + (1  - !>*/* 
far  i ■ 1,2,  ...,  n+1 
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we  determine  8^  to  that 


is  a minimum  which  yields 

(2.3)  h " -(A+1  - ^ 1>/2  + n^X  * 

For  i ■ 1,  2,  ...t  n,  let"*" 

(2JJ  s*  - (S±  ♦ Si+1)/2  . 

From  the  plot  of  the  r vs  s function  the  values  A p ^2*  •••*  ^n 
slf  s2,  ...,  an  (the  midarcpoints  between  the  xeros  of  6)  can  be 

measured. ^ 

Let  Oj.  for  i,j,  - 1,  2,  ...»  n be  such  a real  number  that 

(2.5)  " V<VUi)3  “P  t"aij(sj  " 1 2i)J 

where  Up  Up  are  the  values  of  u for  Sp  s2,  •••>  Bn«  ^ 

equation  (1.21)  were  satisfied  for  s » Sp  r - for  i • 1,  2,  ...,  n, 
then  would  be  the  same  for  i,3  “ 2,  ...»  n.  Since  the  data  con- 

tain errors,  (2.5)  is  solved  for 

(2.6)  a±i  - -[log(rJ/ri)  - (3/2)log(uj/ui)]/(s;J  - s±) 

and  the  least  square  value  a is  given  by 

2 n n 

(2.7)  « ■ rr — rr  51  21  #<i  • 

n(n  - 1)  ^i  i-j+1  1) 

1 If  the  data  are  not  reliable  on  all  channels,  it  may  be  desirable  to 
use  the  above  procedure  with  one  channel,  say  ag,  in  place  of  r for 

the  determination  of  X/2  and  Sp  s2,  ...,  sn,  provided  that  the  spin™ 

of  the  boob  introduces  no  new  xeros  in  ag. 

2 An  approximate  method  of  malrlng  a aero-level  correction  on  individual 
channels  is  given  in  the  appendix. 
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Using  (1.21)  the  least  square  value  for  C is  given  by 


C - ^ X exp  (osj)  • 


Fitted  values  r 


1»  r2»  •••» 


of  the  resultant  acceleration  due  to  yarn 


can  then  be  obtained  by  substituting  into  (1.21)^ 


(2.9) 


rl  " exp  (-aai)  . 


For  each  i,  i - 1,  2,  ...,  n,  use  r^  (1.20),  (1.18),  (1.22)  to  obtain 
a corresponding  value  of  6,  6 and  Kg.  A least  square  value  of  Kg  can 
then  be  obtained  by  averaging  the  values  of  Kg  for  i ■ 1,  2,  n. 

Nov  use  (1.23)  and  (1.21;)  to  obtain  K^  and  Kg. 

Spin.  Using  (1*27)  ve  can  compute  the  Cs  vs  t and  Sn  vs  t functions. 
Equations  (1.28)  ensure  that  Cs  and  Sn  are  estimates  of  cos  y and 
sin*/''  . The  variances  of  Cs  and  Sn  become  large  at  zeros  of  a.  Draw 
smooth  curves  through  the  plots  Cs  vs  t and  Sn  vs  t taking  into  account 
the  above  property  of  the  variances.  (See  Figure  5>.)  Henceforth  ve 
use  these  smoothed  values  of  Cs  and  Sn.  Equations  (1.29)  provide  us 
with  two  estimates  of  y . 

Let  Vc,  Vs  be  the  variances  of  y as  estimated  by  arccos  Cs, 
arcsin  Sn  respectively.  If  these  estimates  were  independent  the  best 
(least  variance,  linear  combination)  estimate  of  y would  be 


(arccos  Cs) 


Jc  ♦ (arcsin  Sn) 

1c  ♦ lA» 


Further  if  Cs  and  Sn  had  equal  variances,  Vc,  Vs  would  be  approximately 

2 2 

proportional  to  l/sin  y , l/cos  y , respectively.  In  any  case  a very 
good  estimate  of  y is  given  by 


(2.10) 


2 2 
Sn  (arccos  Cs)  » Cs_ (arcsin  Sn) 

Sn  + Cs2 


The  rate  of  spin  can  new  be  obtained  by  differentiating  numerically. 


3.  S3  FORMULA  AND  D3FFBBBNTI4L  KJOATKH  FCR  S3  IAW 

The  differential  equation  for  the  plane  yawing  of  a bomb  is  derived 
by  MeShane  as  equation  4.6  of  £lj.  Let  9 be  the  elope  of  the  trajectory 
in  a vertical  plane  and  g be  the  acceleration  of  gravity.  Then  the  abeve 
equation  is 

(3.1)  6 ♦ (pd^/mk2)  (Kjj  ♦ k2K^)6  + (pdKjjuVmk2  )6 

♦ (pd2g/nk2)(XH  ♦ k2Kp)cos  0 - (2g2^2)»in  0 cos  0 • 0. 

(banging  to  arclength  derivatives  and  using  (1.9)  and  (1.10)  this  equa- 
tion becomes 

(3.2)  6"  ♦ (2a  + uVu)6«  ♦ p26  + (ag/u2 )cos  0 

- (2g2/u^)sin  9 cos  0 • 0 . 


Using  the  notation  of  MeShane  we  let  2Q.,  Qp,  Q.  be  the  coefficients  ef 
6',  6,  1 of  (3.2)  and  write  (3.2)  as  ** 


(3.3) 


6"  ♦ 2^6*  ♦ 0^6  ♦ Qj  - 0 


The  complete  solution  of  (3.3)  as  given  by  (7.1)  of  111  is  (where 
a is  a bound  variable  replacing  s and  c^,  c2  are  complex  constants) 

(3*4)  8(a)  * 6^(s)  ♦ ^ exp  [/  V«>d0  + 1 9 (*)] 


* 

°2  exP  §2  Wfo  - i ^P(»)J 


subject  to  the  conditions  that 

(3.5) 


r -vV 


v- v 


(3.6)  f"-0. 

The  second  and  third  terms  of  (3*4)  are  solutions  of  the  homogeneous 
equation  corresponding  to  (3.3)  and  the  first  tern  is  a non-oaclUatory 
particular  solution  of  (3.3).  The  latter  term  is  called  the  yaw  ef 
repose  and  is  given  by  (£.1)  of  [ijj . 
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(3*7)  8^  - (2g2/P^^)«la  8 cos  8 - (flg/p2u2)oo»  • ■ 0, 

For  the  beginning  of  the  trajectory  for  the  T55  and  T$6  bombs  6^  did 

not  exceed  $ x 10"3  radians?  In  general  we  will  make  the 

/ ' 

Assumption  7.  The  yaw  of  repose  can  be  neglected  in  (3.U)  and  (3*2)* 
With  this  assumption  (3*2)  becomes 

(1.7)  6"  ♦ 2 (a  + u,/2u)6l  + p26  - 0 


Substituting  in  the  values  of  Qp  Qg  into  (3*5)  we  find 
(3.8)  <p  1 ■ VP2  “ ^ " su’/u  * (u'/2u)2  - u"/2u  , 

Now 


u'  - u/u  , 

u"  - u/u2  - u2/u3  . 


For  the  T55  and  T£6  bombs  the  following  approximate  values  were  obtained  i 

u - 500  ft. /sec.  a2  - 10‘7/*t.2 

u - 12  ft. /sec.2  P2  - 10"Vft.2 

u - 1.5  ft  ./sec3  nu'/u-lix  10"8/ft.2 

u»  - 2 x 10“2/sec.  (u'/2u)2  - U x 10”10/ft.2 

u"  - 6 x UO^/ft.  sec.  u"/2u  - 6 x 10~9/ft.2 

In  general  we  can  make  the  following 

Assumption  8.  t? t ou’/u,  (u'/2u)2,  u"/2u  can  be  neglected  in  comparison 
with  p2  in  (3.8). 

Using  this  assumption  (3.5)  becomes 


(3*9)  <p * - P 

and  (3.6)  is  satisfied.  Integrating  (3.9)  we  obtain 
(3.10)  f (s)  - ps  + y 


* For  the  values  of  the  ingredients  in  this  estimate  see  the  discussion  of 
Assumption  8 in  Section  3* 
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where  ? ie  a real  constant.  Farther 


(3.11) 


We  can  take  cor  reference  fraae  so  that  the  yaw  6 Is  real*  Then  c^,  Cg 

of  (3.U)  are  c copier  conjugates.  Using  (3.10),  (3.11)  and  assumption 
7 equation  (3.U)  becomes 

(1.8)  6 » ,yu0/u  e“a*(cei^*  ♦ ce-i^*) 

where  c and  % are  complex  conjugates  and  c ■ c1  exp  if,  Z - CgexpC-If). 


U.  RESULTS  AMD  CONCLUSIONS 

As  remarked  in  the  Introduction,  the  experiment  was  not  altogether 
successful,  for  the  aerodynamic  coefficients  obtained  by  the  methods 
discussed  in  Sections  1 and  2 are  for  the  most  part  of  doubtful  accuracy. 

Two  reasons  for  the  poor  results  can  be  attributed  directly  to  the 
electrical  equipment  used.  First,  it  was  necessary  to  measure  drag  from 
accelerometers  of  greater  range  (and  consequently  of  less  accuracy  at 
lower  speeds)  than  had  been  planned.  Second,  the  transverse  accelero- 
meters showed  rapid  and  spectacular  fluctuations  (cf.  Figures  1 and  2), 
which  made  the  analysis  difficult  and  introduced  large  relative  errors 
into  the  computations. 

Because  the  yaws  were  small,  and  damped  out  fairly  quickly,  the 
aerodynamic  coefficients  (except  E^)  and  the  spin  were  composed  only  for 

the  first  part  of  the  trajectory,  at  a Mach  number  of  about  0.6. 

The  following  discussion  of  the  data  for  each  drop  is  given  to 
illustrate  the  difficulties  in  analysis  mentioned  in  the  preceding  para- 
graph, how  some  of  the  difficulties  were  overcome,  and  their  effect 
upon  subsequent  computations. 

Drop  1.  The  accelerometer  corresponding  to  Channel  1 did  not  function 
at  all,  so  the  most  important  data  were  missing.  No  attempt  was  made  to 
reduce  the  data  from  other  channels. 

Prop  3.  The  data  must  be  considered  as  no  better  than  fair.  Although 
Channel  6 displayed  fairly  well  defined  extrema  and  zeros,  Channel  7 
(see  Figure  2)  fluctuated  up  and  down  widly,  rendering  the  data  useless. 
Instead,  Channels  k and  $ were  used.  These  channels  gave  much  better 
curves,  but  the  small  amplitudes  decreased  the  relative  accuracy  of  the 
readings.  Channels  2 and  3 presented  fairly  good  curves. 
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Drop  U.  Because  Channel  7 did  not  function,  the  data  are  useless  fer 
captations  ether  than  Kg,  Kg  and  X.  Channel  6 was  fair}  frsn  it  alone 

X was  determined  and  Kg  was  computed,  in  attempt  to  use  Channels  i*  and 

$ in  place  of  Channels  6 and  7 produced  meaningless  results. 

Drop  2.  On  the  whole,  the  data  are  very  good.  Channel  7 has  sene  wild 
^actuations  downward,  hut  if  these  fluctuations  are  considered  as  error 
and  Ignored,  the  curve  has  well  defined  extrema  and  zeros.  The  other 
channels  are  very  good. 

Drop  5.  The  data  are  only  fairly  good.  Channel  6 is  good,  but  Channel 
7 has  very  wild  downward  fluctuations  - more  than  Channel  7 of  Drop  2. 

The  upper  envelope  of  Channel  7 seemed  to  be  in  phase  with  Channel  6 and 
this  envelope  was  used  to  determine  r.  Use  of  this  envelope,  however, 
makes  it  Inpossible  to  estimate  the  error  in  Channel  7.  Thus  there  are 
unknown  errors  in  Kg,  Kg,  and  Kg  introduced  through  unknown  errors  in  r 

and  a*  Channel  3 is  very  good.  Channel  2 has  wild  downward  fluctuations, 
but  again  the  upper  envelope  was  used. 

Figures  3 and  U show  the  computed  values  of  Kg  vs  Mach  number  for 

the  T££  and  T£6  bombs  respectively.  The  solid  curve  of  each  figure  la 
one  probable  Kg  curve  defined  by  the  data,  faired  in  by  eye.  The 

broken  curves  of  Figures  3 and  U represent  the  error  in  Kg  due  to  an 

error  of  0.02  units  of  gravity  which  is  the  probable  maximum  error  in 
the  smoothed  Sg  vs  t curves.  It  was  assumed  that  the  sere  level,  after 

correction,  is  within  O.Olg  of  the  true  zero  and  that  the  curve  drawn 
represents  the  acceleration  (uncorrected  for  zero  level)  within  O.Olg. 

An  appreciable  error  in  Kg  could  result  from  an  error  in  p.  In  seme 

cases  the  meteorological  data  were  not  taken  until  several  hours  after 
the  drop.  Thus,  no  satisfactory  estimate  of  the  error  in  p can  be  made. 
If  the  determination  of  u were  in  error  by  10  ft./sec.,  at  £00  ft./sec. 
this  would  introduce  into  Kg  an  error  of  only  l&,  a small  fraction  of 

the  error  introduced  by  an  error  in  a^.  As  us  increases,  this  error 

becomes  rapidly  less  Important,  and  hence  Eg  is  relatively  unaffected 
by  errors  in  u. 

Using  the  curves  of  Kg  in  Figures  3 and  U,  Mr.  E.  S.  Martin,  of 

these  laboratories,  has  kindly  computed  the  trajectories  of  a number  of 
bombs  which  were  dropped  in  the  range  bombing  work  at  Edvards  AFB.  Thus 
the  range  and  time  of  flight  could  be  compared  with  the  observed  values. 

It  appears  that  the  values  near  a Mach  number  of  1.1  are  prebably  too 
lav  for  both  the  3,000  lb.  and  10,000  lb.  T£6  bombs. 
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The  computed  values  of  Kg,  Kg,  K^  and  Kg  for  each  drop  art  given  la 

Tibia  1.  Although  It  vaa  lapoaaiblo  to  oatlaat*  the  orrora  la  these 
coefficients,  a few  general  raaarka  oaa  be  made.  The  accuracy  of  Kg  la 

poor  bacanaa  tj  aad  a^  are  snail.  Sven  though  la  general  these  coeffi- 
cients as  given  In  Table  1 are  no  better  than  an  indication  of  the  order 
of  magnitude,  these  for  Drop  2 are  probably  of  useful  accuracy. 

• 

The  angle  of  rotation y and  the  rate  of  ■pin'#'  for  Drops  2 and  3 
are  shown  In  Figures  6 through  9,  the  optical  and  Doppler  values  of  Y 
(adjusted  for  phase  differences)  being  shown  for  comparison.  The  baib 
of  Drop  5 displayed  no  discernible  rotation.  The  spin  rate  8 seconds 
after  release  is  given  in  Table  1. 


Table  1 


Drop  No. 

3 

1* 

2 

5 

Bomb  Type 

T55 

T55 

T56 

T56 

Serial  No. 

6161* 

6169 

6116 

6132 

Date  Dropped 

U Mar  52 

l*  Mar  52 

25  Feb  52 

5 Mar  52 

m (lbs.) 

2965.5 

2989.5 

10932 

10900 

d (ft.) 

2 

2 

3.833 

3.833 

B (lb .-ft.2) 

23083 

22229 

99097 

1019U7 

h (ft.) 

1.55 

5.57 

5.57 

a 

.000296 

1 

.000229 

.000338 

P 

.OOU99 

1 

.001*16 

.00397 

X/2  (ft.) 

630 

7UU 

756 

791 

L (Mach  .6) 

.O8I45 

.081*5 

.01*95 

.01*95 

*M 

2.70 

1.72 

1.23 

1.11* 

*N 

0.71 

1 

1.51 

3.20 

“l 

0.63 

l 

1.1*6 

3.15 

30.9 

l 

7.6 

10.7 

Rafter  8 see.) 
(rad. /sec.) 

1.11 

1 

.175 

0 

Channel  7 did  not  function. 


t* 
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Comparison  with  Cornell  data.  For  purposes  of  comparison,  end 

were  computed  from  Cornell  Aeronautical  Laboratory  [2].  The  Cornell 
report  gars  graphs  of  C^  and  Cg  rs  a (here  a is  the  angle  of  attack  In 

degrees)  far  the  standard  3,000  and  10,000  lb.  bombs.  Now  C and  K are 
related  by 


(U.D 


( \ - - (180/8 ) (dCj,j/dc)0 

[ Kjj  - (180/8)  (dCjj/d«)0 


where  the  subscript  zero  Indicates  that  the  derivative  is  evaluated  far 
zero  degrees  angle  of  attack.  If  and  Cg  are  defined  as  least- 

squares  cubic  equations  in  a,  the  values  of  and  Kg  can  be  computed 

directly  from  (U.l).  The  positions  of  the  centers  of  gravity  of  the 
bombs  dropped  at  White  Sands  differed  by  about  0.1  calibers  from  those 
of  the  standard  bombs  for  which  the  Cornell  data  were  computed.  Hence, 
equations  (U.l)  oust  be  multiplied  by  the  ratio  of  the  c.g.  distances. 
Table  2 is  a comparison  of  our  results  with  the  adjusted  Cornell  results 
for  Mach  number  0.6.  Hie  superscripts  0 and  c indicate  our  results  and 
the  adjusted  Cornell  results  respectively. 


Table  2 


Drop  No. 

3 

U 

2 

5 

V 

2.70 

1.72 

1.23 

1.1U 

V 

2.U0 

2.U2 

1.25 

1.2U 

V 

0.71 

— 

1.51 

3.20 

V 

1.17 

1.17 

1.52 

1.52 

Fran  the  agreement  between  the  results  of  Drop  2 and  those  of  the  Cornell 
data  it  can  not  be  concluded  that  the  scale  effect  is  negligible  for 
this  bomb,  even  though  the  accelerometer  data  are  quite  smooth,  because 
the  sample  of  good  data  is  too  small. 

Conclusions.  -Our  remits  indicate  that  this  method  of  getting  aerody- 
namic  data  is  accurate  and  practical  if  the  following  improvements  can 
be  made: 

(a)  The  electrical  system  can  be  made  to  give  a smooth  record  of 
the  accelerations,  with  fairly  well-known  limits  of  error.  This  is  the 
moat  important. 
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(b)  The  dreg  meaaurementa  can  be  nede  with  two  or  more  accelero- 
meter* of  different  range* , to  obtain  email  relative  error*  with  lew 
drag. 

(o)  Accurate  meteorological  data  can  be  given  near  the  time  of  the 

drop. 

If  I*.  Ig,  .to.,  ar.  wonted  for  Mach  tartar.  orach  hl*ar  than  that 

of  release,  it  may  be  necessary  to  excite  oscillations  late  in  the  flight 
by  eome  device. 


J.  OONLAN 


A.  S.  GALBRAITH 


J.  T.  LEVIS 


AFFBNDHs  AFFRQXIMATB  ZKR0-LE7BL  COERECTION 


If  the  airplane  is  not  accelerating,  the  portion  of  the  longitadinal 
acceleration  curve  before  the  time  of  release  should  be  zero.  With  this 
assumption,  the  smoothed  a^  vs  t curve  to  the  left  of  the  time  of  release 

indicates  the  true  zero  level  of  a^,  and  any  needed  zero  level  correction 
can  be  made. 


If  the  indicated  zero  for  a particular  transverse  acceleration  j is 
in  error,  a corrected  zero  can  be  found  approximately  assuming  the  com- 


ponent Sj  behaves  like  r. 


Let  a^s^ 


and 


s^)  be  the  corrected  and 


indicated  accelerations,  respectively,  at  8^.  We  attempt  to  determine 
kj  so  that 


(A!)  aj(®i)  " *V8i)  “ kJ  * 


Using  (1.18)  and  the  equal  spacing  of  s^,  s2,  ...,  sn  we  have 

<*2)  (V“M]‘V2 

■ [vi>i.3'3A  ‘iW'jW 

for  i ■ 1,  2,  n - 2.  Making  the  approximation 


uiui+2 


in  (A2)  we  obtain 

(A3) 

or  using  (Al)  in  (A3), 


for  i ■ 1,  2,  n - 2.  The  least  square  solution  of  (AU)  for  kj  is 

n-2  n-2  ~ 

(W)  k3  ’ j£  piV'1^  pi 

where 

pl  " - »j(«i)*j(«l.2> 

and 

^i  “ 2a^Bi+l*  ’ a^si*  " a^ai+2^* 

2U 


Liat  of  Qpfcols 


V *2 


acceleration  as  recorded  an  channel  j,  j ■ 1,  1,  ...,  7. 
trance  acceleration,  a2  - a62  ♦ a/, 
transverse  «c«ent  of  inertia  of  bcab. 


, sign  corresponding  to  that  of  a. 


l(/^re 


dimeter  of  bcab. 


acceleration  due  to  gravity. 


distance  from  center  of  gravity  of  bonb  to  acceleroneters. 
radius  of  gyration  of  bonb 
mass  of  boob. 

the  mraber  of  Intervals  defined  by  successive  zeros  of  6. 

Vi  2 2 / 2 2 

a6  + *7  “ V *2  + *3  ' 

distance  along  the  trajectory. 

“*7/±'1^6r  ♦ i if , sign  corresponding  to  that  of  a. 
tine. 

air  speed  of  bcab. 

(pd2/2ak2)(k2KL  + Kh). 

*Y  pd 3km/b. 


two  coaplex  terns  of  6 in  equation  (1*8). 
slope  of  the  trajectory  in  a vertical  plane, 
distance  between  successive  zeros  of  6. 
air  density. 

$a  ♦ (real  constant), 
argment  of  6^,  precession. 
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*P2  ■ argument  of  S>2,  nutation. 

ty/  • angle  of  axial  rotation  of  the  bomb. 

<o  ■ angular  velocity  of  the  longitudinal  axia  of  the  bomb. 

A dot  (*)  denotes  a derivative  with  respect  to  time,  and  a prime  (') 
denotes  a derivative  with  respect  to  arc  length. 

AERODYNAMIC  COEFFICIENTS 


Kp  » drag  coefficient. 

Kg  - damping  moment  coefficient. 

■ lift  coefficient. 

Kg  - restoring  moment  coefficient. 
Kg  - normal  force  coefficient. 


AERODYNAMIC  FORCES 

2 2 

D - drag  force  ■ Kppd  u . 

H * damping  moment  « -KgpdKico. 

L - lift  force  - KLpd2 3u26. 

M ■ restoring  moment  ■ -Kgpd^u26. 

2 2 

N ■ normal  force  • K.Tpd  u 


The  following  special  symbolism  regarding  arc  length  and  accelera- 
tion is  used  in  Section  2. 

1.  Capital  letters  refer  to  values  of  the  functions  at  the  zeros 
of  8. 

2.  Small  letters  refer  to  values  of  the  functions  at  the  extrema 
of  6. 


3.  Script  letters  refer  to  observed  or  unfitted  values. 
1*.  Printed  letters  refer  to  fitted  values. 
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